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Using the luminosity function and the apparent GRB rate inferred from the spectral 
peak energy-luminosity relation, we investigate the absolute GRB formation rate, taking into 
account the effects of both the jet-luminosity evolution and the jet opening angle evolution. 
In the case that there is no jet opening angle evolution, the jet-corrected luminosity for 
high-redshift GRBs is much larger than typical values for low- z (z ~ 1) bursts. On the other 
hand, if the jet-luminosity does not evolve with time, the jet opening angle for high- z bursts 
should be much smaller than that for low-z GRBs. Therefore, it is preferable to take into 
account both evolution effects. We also estimate the local GRB event rate in a galaxy of 
^ . ~ 10- 7 -10" 5 yr" 1 . 

81. Introduction 

The luminosity function and the formation rate history of gamma-ray bursts 
(GRBs) are important for understanding the nature of these events. The formation 
history of GRBs provides constraints on the population of progenitors. If we know 
the absolute value of the GRB rate, we can study their relation to other objects 
potentially related to GRBs. It is fairly certain that long GRBs arise from relativistic 
jets with typical opening half-angles of ~ 0.1 rad.E3B3 For this reason, determining 
the corrections due to jet collimation and/or jet structure is essential in the study 
of the GRB rate. The luminosity function may directly reflect the properties of the 
progenitor, such as the jet structure, the jet kinematics, and so on. 

GRBs could be useful tools to investigate high redshift objects. The connection 
between long duration GRBs and supernovae is strongly suggested, and it is known 
that at least some GRBs arise from the collapse of massive starsl^'^U^U^ The 
universe began to reionize at z ~ 17p^ which implies that the first generation objects 
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had already been formed by at least z ~ 20. Because they are massive, they imme- 
diately exploded and produced GRBs at a redshift of z ~ 20. Indeed, some of their 
prompt emissions and afterglows are bright enough to be actually detected The 
absolute rate of prompt GRB emissions may provide some information concerning 
the cosmic star formation history I22H22I Also, the afterglows of high-z GRBs may 
provide information concerning the cosmic reionization history ED03E3JniJ3 

Many authors have attempted to derive the luminosity function and absolute rate 
of GRBs. Schmidt 341 derived the luminosity function and the absolute value of the 
local GRB rate so that (V/V maiX ) and the observed event rate of BATSE bursts are re- 
produced. Later, he obtained the (U/V max )-hardness relation for BATSE GRBs that 
was used to determine the luminosity function and the absolute local GRB rateP^ 
The value he derived is ~ 0.5 Gpc _3 yr~ 1 . However, in that derivation, the geomet- 
rical corrections were not considered. The luminosity function and the formation 
rate history of long GRBs have been simultaneously derived using some luminos- 
ity indicators, such as the luminosity- variability relation^ ^3 and the luminosity-lag 
relationP^J^ Recently, it was found that the spectral peak energy is correlated 
with the peak luminosity with actually known distances of up to z ~ 4.5.^ The 
correlation is tightest in such kinds of luminosity indicators Using the peak 

energy-luminosity relation as a luminosity indicator, Yonetoku et al.^ investigated 
the apparent GRB luminosity function (see § 2.1) and the GRB formation rate for 
redshifts z < 12 and suggested that the GRB formation rate rapidly increases up to 
z ~ 1 and then continues to increase toward higher redshifts z ~ 12. This behavior 
of the formation history is somewhat similar to those derived using other luminosity 
indicators EHsaJaa 

There is evidence suggesting that the (isotropic-equivalent) luminosity evolution 
of GRBs takes the form oc (1 + z) 2,6 P^ which is quite similar to that of quasars .^^^ 
The luminosity evolution for GRBs has been independently reportedP^''^)^'^ 
The origin of the luminosity evolution is not yet clear. The jet-corrected luminosity, 
Lj = Lfb, where = 0j 2 /2 is the beaming factor and 9j is the jet opening half- 
angle, clusters very near the standard valued The same is true for the jet-corrected 
energy da These facts suggest that the luminosity evolution may result from 
either jet-opening angle evolution^ or jet-luminosity evolution. 

The effect of jet collimation is investigated explicitly in Ref. I14|) . where the 
luminosity function is studied assuming that its functional form is a double power- 
law and that the jet-corrected luminosity Lj is constant. For the structured jet case, 
general formulae to calculate the luminosity function and/or the local GRB rate 
from the angular distribution of kinetic energies in the jet has been presented .^J^ 
One crucial assumption in these works is that the apparent GRB rate is proportional 
to the cosmic star formation rate. Furthermore, in Ref. I14[). luminosity evolution is 
not considered. 

In this paper, using the luminosity function and the apparent GRB rate derived 
in Ref. I41|) . we study the absolute GRB formation rate, taking into account the 
effects of both the jet-luminosity evolution and the jet opening angle evolution. A 
general formula to calculate the true GRB rate is presented. The most important 
point in the derivation is the geometrical correction of the jet opening angle, 9j. 
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Instead of using the distribution function for 9j, we consider an approach similar 
to that used in Ref. I14|) . The fact that the jet-corrected luminosity Lj clusters in 
narrow ranee enables us to calculate the jet-correction factors without directly 
taking into account the opening angle distribution. We also study the opening angle 
evolution. If we assume that the origin of the luminosity evolution is the jet opening 
angle evolution, the high-redshift GRBs are collimated much more tightly than the 
nearby events. Throughout the paper, we assume a flat, isotropic universe with 
Q M = 0.32, Q A = 0.68, and h = 0.72. 

§2. The absolute GRB formation rate 

2.1. Apparent luminosity function of GRBs 

The apparent luminosity function of GRBs, that is the rate of actually detected 
GRBs per unit comoving volume at a redshift z and a peak luminosity L, is written 

<P(L,z)=p(z)4>(L/g(z))/g(z) , (2-1) 

where p(z), g(z), and 4>(L') are the GRB formation rate, the luminosity evolution, 
and the local luminosity function, respectively!*'*^'''^' In this paper, we adopt the 
functional forms of these quantities given below, as derived in Ref. 41). 
The GRB formation rate is given by 

< z < 1 , 

1 < z < z max , (2-2) 

2-max Z ; 

where po = p( z = 1) is the normalization constant, and z max is the epoch at which the 
first GRBs formed. Here we assume a broken power-law form. However, the results 
obtained in the present paper are almost unchanged if other smooth functions are 
adopted. 

The luminosity evolution is assumed to take the form 

9^ = {^) a . (2-3) 

where, throughout the paper, we adopt a = 2.6, as derived in Ref. I4*T|) . 

The local luminosity function 4>(L'), where V = L/g{z) represents the peak lu- 
minosity after removing the luminosity evolution effect, is derived from the universal 
cumulative luminosity function ip(L') as 

<KIJ) = ~W) , (2-4) 
and we assume the following form of ip(L') in the range L' min < L' < L' mSLX : 

W) = (L'/K) A [1 + (L'/L:) s ] {B ~ A)/s . (2-5) 

Here A and B are the power law indices in the asymptotic forms, and s describes 
the smoothness of the transitions between the high and low L' regimes. We set 



p{z) = Po 



[(1 + z)/2] 

[(1 + *)/2] 
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these parameters as A = —0.29, B = —1.02, and s = 2.6, in order to reproduce 
the results presented in Fig. 6 of Ref. I41[) . Therefore, the break point of ip{L') is 
L'^ = 1.1 x2"x 10 51 erg s _1 , because the functional form of g(z) adopted in this 
paper has an extra factor of 2 a compared with that used in Ref. 14 1|) . Finally, the 
normalization constant ipQ is chosen so as to satisfy 

nax 

4){L')dL' = 1 . (2-6) 

n 

As suggested in Ref. 4TJ), we assume the universality of the local luminosity 
function <p(L'). Then, both L' max and L' min should be independent of z; that is, we 
can write 

Ljoaxiz) = L' max x g(z) , f2-71 
Lmm(z) = L' min x g(z) . 

Hence both the maximum and the minimum values of the luminosity evolve with g(z). 
In the fiducial case, the values L' max = 1 x 10 54 erg s _1 and L' min = 1 x 10 49 erg s" 1 
are adopted.^ 

We determine the normalization constant po as follows. The actual event rate 
of BATSE-observed bursts, with observed peak fluxes larger than -F b s > is given by 

N ohs {> F ohs ) = Cexp / dz-—— dL<P(L,z) , (2-8) 

J l + zdz J LlUz) 

where V, L\[ m , and zn m are the comoving volume, the limiting flux for given z, and 
the upper bound on the redshift, respectively. The mean exposure efficiency Cexp is 
approximately 48% for the BATSE experiment. 301 The limiting flux can be written 
as 

L hm {z) = A7rd L {z) 2 x F obs , (2-9) 

where di{z) is the luminosity distance. In Ref. I41|) . 689 BATSE bursts were con- 
sidered. Their observed peak fluxes in the 30-10000 keV band are larger than 
2 x 10~ 7 erg cnr 2 s _1 and their redshifts are smaller than 12(= Zu m ). The cor- 
responding observation period is 9.0 yrs, because they selected the GRB sources 
from ID# 105(1991.111) to ID# 8121(2000.147) (see Table 2 of Ref. EJ)). Hence 
we find iV bs(> 2 x 10~ 7 erg cm _2 s _1 ) ~ 76.6 events yr _1 . We can determine the 
value of the normalization constant po needed to reproduce this value. The results 
are shown in Table [I] It is found that po changes by a few percent if L' max ranges 
between 1 x 10 53 and 1 x 10 55 erg s _1 . The dependence of L' min is estimated as 
follows. Equation (|2-6j) . together with Eqs. (|2-4|) and (|2-5|) . gives the approximate 
form Vo ~ (L' m jL'J- A . Then, Eq. reads p cc Vo" 1 oc (L' m jL^) A . Therefore, 
as shown in Table [JJ when L' min becomes 10 times larger (smaller), po becomes about 
2 times smaller (larger). Finally, we find that p$ changes by a factor of 2 when the 
parameters in the luminosity function are varied in ranges corresponding to 90% of 
the uncertainties on the observational results derived in Ref. I41j) . 
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Table I. The values of po = p(z — 1) (in units of Gpc 3 yr for various luminosity ranges [see 
Eq. fZty ]. 





1 x 10 53 


rl a) 
^max 

1 x 10 54 


1 x 10 55 




1 x 10 4S 


1.69 


1.63 


1.62 


rl a) 

■"mm 


1 x 10 49 


0.86 


0.84 


0.83 




1 x 10 50 


0.44 


0.43 


0.43 



a) The maximum and minimum luminosities in units of erg s 1 [see Eq. (12-71 1 . 



We obtain p(0) ~ 0.01 Gpc _3 yr _1 and p(l) ~ 1.0 Gpc _3 yr _1 in this paper. 
Let us compare these with previous results. From the results obtained in Ref. IHJ, 
in which the functional form of p{z) is determined so as to reproduce both the 
observed peak flux and redshift distributions, the values p(0) ~ 0.04-0.4 Gpc _3 yr _1 
and p(l) ~ 0.2-3.8 Gpc _3 yr~ 1 can be inferred. We thus find that our results are 
roughly consistent with those derived in Ref. EJ)- However, SchmidtP^ estimated the 
local GRB rate density to be p(0) ~ 0.5 Gpc _3 yr _1 , which is more than one order 
of magnitude higher than our value. This difference may result from an insufficient 
amount of data and low resolution in the redshift range < z < 1. The form of p(z) 
we employ might be steeper than the actual one. Therefore, because the function 
p(z) is normalized at z = 1 in this paper, the local GRB rate that we derived at 
z = might be smaller than the actual value. 

2.2. GRB formation rate corrected by the jet opening angle 

We consider the true rate of GRBs which actually occurred in the redshift range 
Z\ < z < Z2- For this purpose, it is necessary to consider the geometrical correction 
of the jet opening angle, 9j. Although the jet structure is still unknown, we employ 
a uniform jet model. Taking into account the correction due to jet collimation, the 
cumulative GRB rate is calculated as 

f Z2 1 dV /' imax ( 2 ) 
N( Zl ,z 2 )= dz-—— dL$(L,z) 
Jzi 1 + a dz J LmUz) 



j dL 3 P{Lj,z)h{L,hr X , (2"10) 



where Lj and fb are the peak luminosity confined in the jet and the beaming factor, 
given by 



fb(L, Lj) = 1 — cos 9j(L, L 



3 1 



= ^ (< 1) , (2-H) 

respectively. The quantity Lj is assumed to be clustered in a narrow range around 
the standard valued and its distribution is approximately log-normal: 



V ZTTcrLj 



1 /. _ . , A 2 ' 
2^2 



log Lj - log Lj (z) 



(2-12) 
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where log Lj (z) and a are the mean and standard deviation of the luminosity distri- 
bution, respectively. In this paper, the evolution of the jet-corrected luminosity is 
assumed to take the simple form 

1 + z^ 



L j {z) = L jl \—-\ , (2-13) 

where j3 is the slope index of the jet-luminosity evolution. We discuss in the 
validity of our assumed functional forms of P(Lj, z) and Lj. 

Substituting Eqs. (J2TJ, (EZHJ) , (l2~T2l . and (t2^|) into Eq. (|2~TU1) . we obtain the 
cumulative GRB rate as 

dZ —zTz p{z) ^* ' (244) 



where 



Lji \ 2 



\ I max 

J y dL'<j>(L')L' (2-15) 



is the luminosity weighted beaming factor, and 



J dL j P(L j ,z)L 



-l 

3 



10 i CT 2lnlO 



(2-16) 



is the mean jet-luminosity. 

We choose the value of (Lj) z=1 = L^/IO^ 2 lnl ° so that (/ 6 -1 ) 2 is 200 at z = 1. 
In this case, the mean jet opening half-angle, defined by (0j) z = [(/b" 1 ) /2] -1 ^ 2 , is 
equal to the typical value of 0.1 rad at z = 1. Then we find 



/r A L ii 
\^/ z =l - 10 | CT 2 lnl0 

= 2.5 x 10 49 erg s" 1 f ^^sl ) (2-17) 

with the fiducial parameters (L^ = 1 x 10 54 erg s _1 and L' min = 1 x 10 49 erg s _1 ). 
The value of {Lj) changes by only a factor of 3 when L' max and L' min vary between 
1 x 10 53 and 1 x 10 55 erg s _1 and between 1 x 10 48 and 1 x 10 50 erg s _1 , respectively. 
Substituting Eq. (|2~T7)l into Eq. (jHHJ), we derive 
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l+z 



Fig. 1. Luminosity- weighted beaming factor 
(fb) as a function of the source redshift 
z. The solid, dashed, and dotted lines cor- 
respond to P = a, a/2, and 0, respectively. 
In the case f3 = a (solid line) , the beaming 
factor is independent of z, with the value 
(ft 1 ) = 200. 



10" 



A_ 
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V 
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P = a 
P = a/2 
p = 



l+z 
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Fig. 2. Mean opening half-angle (9j) as a 
function of the source redshift z. The solid, 
dashed, and dotted curves correspond to 
P = a, a/2, and 0, respectively. In the 
case p = a (solid line), the mean jet open- 
ing angle is independent of z, with the value 
(9j) = 0.1 rad. 



and 

W. = ai(i±i) 

where 7 = (/3 — a)/2. 

Figures ^ and [21 display the 
luminosity- weighted beaming factor, 
(f^ 1 ) , and the mean opening half- 
angle, (9j) z , as functions of z. The 
value of (3 is uncertain. As seen in Ta- 
ble [nj we consider three cases: f3 = a 
[no jet-opening angle evolution), (3 = 
a/2 [intermediate), and (3 = {no jet- 
luminosity evolution). In the case (3 = 0, 
we find </ 6 _1 ) 2 ~ 3 x 10 4 and (0 j ) z ~ 
9 x for z ~ 12 in the fiducial case. The latter is about 10 times smaller than 
the typical value for z ~ 1. By contrast, when (3 = a = 2.6, the beaming factor 
is independent of z, and we have (f^ 1 ) = 200. However, the mean jet-corrected 
luminosity Lj increases with z [see Eq. (I2~T31) ] and 4(12) ~ 1.3 x 10%(1), which 
is also unnatural. Hence, intermediate values of (3, in the neighborhood of (3 ~ a/2, 
for which both Lj and (f^ 1 ) depend on z, may be preferable. 

Figure |21 plots the cumulative GRB rate N(0,z) as a function of z, which is the 
GRB frequency at redshifts smaller than z. When we vary (3 from a to 0, -/V(0, z max ) 
increases by approximately one order of magnitude, and the deviation becomes larger 
for higher redshift. 



Table II. Models. 





Evolution 


a 


Jet-luminosity evolution 


(no jet- opening angle evolution) 


a/2 


Inter- mediate 





Jet-opening angle evolution 


(no jet-luminosity evolution) 
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Fig. 3. The cumulative GRB rate N(0, z) [see 
Eq. I|2-14|l ]. The solid, dashed, and dotted 
curves correspond to /3 — a, a/2, and 0, re- 
spectively. These are frequencies of GRBs 
that actually occurred at redshifts from 
to z. 




1 10 
1+z 



Fig. 4. The comoving GRB rate density n(z). 
The solid, dashed, and dotted curves cor- 
respond to j3 = a, a/2, and 0, respectively. 
In the case /3 = a, n(z) is simply propor- 
tional to p(z). 



Finally, the true comoving GRB rate density is given by 

= P(*){/|T 1 } = ■ (2-20) 

Figure |I] presents the results. In the case j3 = a, the value of {f^ l ) z does not depend 
on z, and thus n(z) is simply proportional to p(z). It is also seen that when the 
jet opening angle evolution is considered (/? < a), the GRB rate at high redshift 
becomes much larger than that at low-z. 



§3. Discussion 



We have estimated the absolute GRB rate taking into account the effects of 
jet collimation. The normalization constant, po, was chosen so as to reproduce 
the BATSE detection rate of bright GRBs with observed fluxes larger than 2 x 
erg cm _2 s . In this paper, we specifically considered for the first time both 
the jet luminosity evolution and the jet opening angle evolution. Both evolution 
effects are important because when only one of them is taken into account, extreme 
results are obtained in the high redshift regime, say, z > 10. When the mean jet 
opening angle does not depend on z {(3 = a; no jet-opening angle evolution), the 
jet-corrected luminosity Lj for high-redshift GRBs is much larger than the typical 
values for z ~ 1 bursts. On the other hand, when Lj is independent of z (j3 = 0; 
no jet-luminosity evolution), the jet opening angle for high-z bursts should be much 
smaller than that for \ow-z GRBs. Therefore, the intermediate case (j3 ~ a/2), in 
which both Lj and the jet opening angle varies with z, is most preferable. It is 
also found that for any value of (3, the true GRB rate n(z) does not decrease for 
high-z as long as z < 12, and might not be simply proportional to the cosmic star 
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1 + z 

Fig. 5. Distribution of the jet-corrected lumi- 
nosity Lj = fbL (with fb = 1 — cos Oj , 
where Oj is the jet opening half angle) 
as a function of the redshift. The red- 
shift is estimated using the spectral peak 
energy-peak luminosity relation (the Yo- 
netoku relation),^ and Oj is estimated^ 
using the Yonetoku relation and the spec- 
tral peak energy-jet corrected energy rela- 
tion (the Ghirlanda relation)?^ The solid 
curve represents the truncation of the up- 
per bound of Lj that is caused by the flux 
limit f| im = lx 10~ 6 erg cm _2 s _1 . 




Index of Jet Luminosity Evolution 



Fig. 6. The best index of the jet luminosity 
evolution measured using the r-statistical 
method. When we assume Lj evolution 
as a function of (1 + z) 13 , (3 = 2.05 is the 
best value, at which the data correction de- 
gree becomes zero. We also plot upper and 
lower bounds of the index corresponding 
to 3<t from the best value. These bounds 
are given by = 2.58 and (3 = 1.30. The 
value (3 = lies approximately 8a from the 
best value, and hence we reject this case in 
which Lj -evolution does not exist. 



formation rate. These facts might imply that the metalicity affects the properties 
of the GRB progenitor. For example, in the high-redshift epoch, because of the low 
metalicity, the opacity of a star was small and the stellar envelope may not be blown 
off completely. Then, it might be difficult for the GRB jets to break out of the stellar 
surface, resulting in small opening angles. 

Recently, the existence of opening angle evolution of the form oc (1+z) 7 with 7 = 
— 0.45lo is was suggested in Ref. 42). Combining this with the result a = 2.6l 2 
derived in Ref. ED, we find (3 = a + 2 7 = 1.70+^;|f , and thus j3 ~ a/2. The jet 
luminosity evolution is seen more clearly in Fig. |SJ Here we use the data set analyzed 
in Ref. 42J). Also, the redshifts and opening angles are estimated in the same manner 
as in Ref. I42|) . In this way, we can obtain the jet-collected luminosity L,- as a function 
of (1 + z). Using the generalized tau statistical methodp^^'^^'E^ we obtain 
(3 = 2.05+0:75 (see Fig. EJ), which is consistent with the value f3 = 1.70t° [H, which is 
predicted using the combination of the luminosity evolution and the jet opening angle 
evolution. Furthermore, Fig. shows that Lj is clustered in a narrow range around 
the mean value for each z and that the distribution of Lj for a given bin of 1 + z 
roughly takes the form of a log-normal with a logarithmic mean a of approximately 
0.6. [When we adopt a ~ 0.6 and </ 6 _1 ) z=1 = 200, we find L jX ~ 3 x 10 49 erg s" 1 
from Eq. (|2T7|) .] We believe that a is independent of z, as assumed in this paper, 
but, the lack of statistics prevents us from carrying out a more detailed study of this 
point. Of course, our method to estimate Lj contains many uncertainties, and hence 
more direct estimations of z and 6j are necessary. If the Swift satellite collects data 
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from many events with spectroscopically measured redshifts and jet opening angles 
determined from the wavelength-independent achromatic break in the afterglow light 
curves, our understanding of evolution effects will be tested. 

Using the local GRB rate, the event rate per galaxy can be estimated as 

rcRB = n(z)/n gsd 

~ 1 x 10- 5 yr- 1 galaxy^ 1 ( - ( " gaI „ ) ' , (3-1) 

y S y VlO 2 yr- 1 Gpc - V VlO 7 Gpc - V ' V ' 

where n ga i is the total number density of spiral and elliptical galaxies, which is 
typically ~ 1 x 10 7 Gpc -3 P^ In the case (3 = a/2, we obtain tgrb = 2 x 
10 -5 yr -1 galaxy -1 for z = 1 and tgrb = 1 x 10 -7 yr -1 galaxy -1 for z = 0. 
Sokolov^S independently estimated the result tgrb ~ 5x 10 -8 (A -1 ) yr -1 galaxy -1 , 
which is roughly consistent with our result. Our values are comparable to the local 
rate of core-collapse (energetic Type Ib/Ic) SNe, ~ 10 -6 -10 -5 yr -1 galaxy -1 ,™ al- 
though the rate for type Ib/Ic SN at high redshift is not known. Even if collimated 
ultra-relativistic outflow and relativistic beaming effect prevent us from observing 
prompt GRB emissions far from the jet axis, the emissions from associated energetic 
core-collapse SNe are nearly isotropic. Hence, the expected SN rate can be estimated 
as 7V(0, 2sn)- The Supernova/ Acceleration Probe*^ will detect SNe with redshifts of 
up to zsn ~ I- 7 , and thus we expect N(0, zsn = 1.7) ~ 3 x 10 4 events yr -1 . 
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